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ABSTRACT 

Using the Expanded Very Large Array, we have mapped the vicinity of the Orion BN/KL Hot Core 
with sub-arcsecond angular resolution in seven metastable inversion transitions of ammonia (NH3): 
(J,K)=(6,6) to (12,12). This emission comes from levels up to 1500 K above the ground state, enabling 
identification of source(s) responsible for heating the region. We used this multi-transition dataset to 
produce images of the rotational/kinetic temperature (Tj-ot/Tkin) and the column density Ncoi of NH3 
for ortho and para species separately and on a position-by-position basis. We find Trot and Ncoi in 
the range 160—490 K and (1—4) x 10^^ cm~^, respectively. Our spatially-resolved images show that 
the highest (column) density and hottest gas is found in a northeast-southwest elongated ridge to the 
southeast of Source I. We have also measured the ortho-para ratio of ammonia, estimated to vary 
in the range 0.9-1.6. Enhancement of ortho with respect to para and the offset of hot NH3 emission 
peaks from known (proto)stellar sources provide evidence that the NH3 molecules have been released 
from dust grains into the gas-phase through the passage of shocks and not by stellar radiation. We 
propose that the combined effect of Source I's proper motion and its low-velocity outfiow impinging 
on a pre-existing dense medium is responsible for the excitation of NH3 and the Orion Hot Core. 
Finally, we found for the first time evidence of a slow (^5 km s~^) and compact (■^lOOO AU) outfiow 
towards IRc7. 

Subject headings: ISM: individual objects (Orion BN/KL) — ISM: molecules — ISM: abundances 



1. INTRODUCTION 

Orion BN/KL is the closest high-rn ass star forming re- 
gion (414±7 pc; iMenten et al.l l20Q7t ) and a compelling 
target for studying how massive young stellar objects 
(YSOs) form and interact with their surroundings. Ob- 
jects of considerable long-term interest in BN/KL are the 
highly embedded radio Source I, believed to be a massive 
YSO (|Reid et al.l 120071: IMatthews et al.ll2010D . and the 
Hot Core >1" away, which is a r i ch so urce of molecular 
emission (e.g.. iGenzel fc Stutzkil 119891 ). In dust emis- 
sion, the compact source SM Al is visible 2"-3" from 
Source I ([Beuther et al.l I2004D. which also shows rich 
molecular emission (jBeuther et alll2005[) . On the other 
hand, only weak molecula r thermal emission h as been 
detected toward Source I fBeu ther et al.l |2005() . How- 
ever, its surrounding mass fiows excite strong H2O and 
SiO maser emission, permitting in vestigations of the 3-D 
gas dynamics at 10-10 00 AU radii (jGreenhill et al.ll2004al: 
IMatthews et al.ll2010D . 

Despite intensive radio and infrared (IR) study, the 
primary heating sources for BN/KL 10^ Lq) are still 
unknown. A strong mid-IR source, IRc2, has been long 
suspected to be the dominant energy source, but when 
observed at sub- arcsecond resolu t ion, i t breaks up into 
multip le peaks (jPougados et al.l 119931 : iGreenhill et al.l 
l2004bD . Moreover, mid-IR observations have recently 
shown that IRc2 is not self-luminous, but is illuminated 
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and heated by Source I (jOkumura et al.|[201lD . While it 
has been proposed that Source I is also (externally) heat- 
ing the Hot Core (jHermsen et al.l[r988( ). direct evidence 
is lacking and the kinetic temperature (Tkin) of molecu- 
lar gas within 1" from Source I is not known. Tackling 
this uncertainty requires studies at radio rather than IR 
wavelengths, in order to penetrate the high column densi- 
ties towards the Hot Core, and measurements of thermal 
rather than maser lines. 

Ammonia (NH3) is well suited to measuring Tkin over 
densities > 10^ cm~^. In particular, metastable transi- 
tions ( J=K) are interesting because they are coUisionally 
excited, thus different (optically-thin) transitions can be 
used to determine Tkin y ia a rotational diagram anal- 
ysis ()Ho fc TownesI 119830 . Inversion transitions from 
(J, i^)=(8,8) to (14,14), with upper state energy levels 
'-~^800-2000K above ground, have been previously de- 
tected toward BN /KL with the Effe lsberg-lOOm tele- 
scope (30" beam; IWilson et al.l 119931 ) , confirming the 
presence of hot molecular gas in the region but not lo- 
calizing it. IWilson et al.l ()2000[ ) mapped the Hot Core 
region in the {J,K) =(4,4) and (10,9) inversion transi- 
tions (^^200 K and ^^1350 K above ground, respectively) 
with >1" resolution. (10,9) is a nonmetastable transition 
of ortho-NH3 and (4,4) is a metastable transition of para- 
NH3 , hence the Boltzmann analysis required assuming an 
ortho-para ratio, preventing an accurate measurement of 

Tkin ■ 

In this Letter, we present new observations of 
(metastable) high-J inversion lines of NH3 in Orion 
BN/KL. Using the Expanded Very Large Array (EVLA; 
IPerlev et al]l20Tll) . we imaged at sub-arcsecond reso- 
lution seven NH3 lines with energy levels high above 
the ground state (equivalent to 400-1600 K), from 
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{J, K)={6fi) to (12,12). The multi-transition measure- 
ments enabled estimation of the temperature, density, 
and velocity field of (hot) molecular gas and hint at a 
possible excitation mechanism of the Hot Core. 

2. OBSERVATIONS AND DATA REDUCTION 

Observations of NH3 were conducted using the EVLA 
of the National Radio Astronomy Observatory (NRAO). 
By using the (new) broadband EVLA K- and Ka-band 
receivers, we observed a total of seven metastable inver- 
sion transitions of NH3, from {J, K) — {Q,6) to (12,12) at 
1.3 cm (25-31 GHz). Transitions from (6,6) to (11,11) 
were observed separately in 3h tracks on several dates 
using the C configuration. The (12,12) dataset was ac- 
quired in the hybrid CnB configuration. Table [1] summa- 
rizes the observations. Each transition was observed us- 
ing a 4 MHz bandwidth (--40 km s"^ at 30 GHz) consist- 
ing of 256 channels with a separation of 15.6 kHz. Typ- 
ical on-source integration was 1.5h. Each transition was 
observed with "fast switching" , where 80s scans on-target 
were alternated with 40s scans of the nearby (1-3°) QSO 
J0541-0541 (measured flux density --0.7 Jy). We de- 
rived absolute flux calibration from observations of 3C 48 
{F^ = 1.1 Jy) or 3C 147 {F^ =1.4-1.6 Jy), depending on 
the epoch, and bandpass calibration from observations 
of 3C 84. 

Using the Astronomical Image Processing System 
(AIPS) task IMAGR, we imaged the BN/KL region with 
cell size 0''1, covering a 50" field. We fitted and removed 
continuum emission from spectral cubes using AIPS' IM- 
LIN. The data were processed both with and without 
velocity-smoothing, resulting in velocity resolutions of 
0.15-0.19 km s~^ and 0.6-0.8 km s~^, respectively (de- 
pending on transition). The following analysis is based 
on the smoothed dataset. To match the angular resolu- 
tion of different configurations, we produced maps set- 
ting the "ROBUST" weighting parameter to TZ—0 for 
the transitions observed with the C-array and to TZ=5 
(natural weighting) for the (12,12) transition observed 
with the CnB-array. We restored all images with a 
beam of size 0''85x0"7 (P.A.=0), approximately the av- 
erage size among different transitions (see Table [T]). Af- 
ter smoothing, the RMS noise per channel was typically 
~1.6 mJy beam~^. The (8,8) transition falls at the band- 
edge of the Ka-band receiver where performance was not 
optimal, resulting in lower-quality data compared with 
the other transitions. Hence, we excluded it from the 
following analysis. Moreover, inspection of data in post- 
processing revealed a 30% offset in the absolute flux scale 
for the (9,9) dataset, possibly due to instability of the 
Ka-band system during the observations. We accounted 
for the systematic error by scaling the flux density by 
30% in the flnal image of the (9,9) transition. 

3. RESULTS 

For the flrst time we have mapped and localized the hot 
NH3 gas from metastable transitions (6,6) up to (12,12) 
with sub-arcsecond resolution towards BN/KL. We ob- 
served three ortho (6,6; 9,9; 12,12) and three para (7,7; 
10,10; 11,11) transitions so as to bracket an upper-state 
energy of 1000 K without exceeding 2000 K. 

For each transition, we produced spectra by mapping 
each spectral channel and summing the flux density in 
each channel map (Figure [l}. Multiple transitions show 
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Figure 1. Spectra of various NH3 inversion transitions observed 
toward Orion BN/KL with the EVLA. The velocity resolution is 
0.7 km s~^ . The fiux density is integrated over the Hot Core region 
and the radial velocities are with respect to the local standard rest 
(LSR). The vertical dashed line indicates a velocity of 6 km s~^. 
Note that different flux density scales are adopted for different 
transitions. 

similar line proflles and central velocities ('^5-6 km s^^), 
while velocity widths, a^, systematically decrease from 
(6,6)^5 km s~^ to (12,12)'^2 km s~^, as determined 
from single-Gaussian fits (Table [1]). By comparing with 
spectral profiles from the Effelsberg-lOOm telescope (as- 
sumi ng an equivalence relation K=0.86 Jv: iWilson et al.l 
Il993f ). most of the single-dish flux density for transitions 
from (9,9) to (12,12) was recovered by the EVLA. This 
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Table 1 

Parameters of Observations. 



Transition ^ 
(J,K) 


t^rcst 

(MHz) 


(K) 


Date 

(yyyy/mm/dd) 


EVLA 
Receiver 


Synthesized Beam 

eui") X e^{") 


RMS 
(mjy/beam) 


^ peak 

(Jy) 


Fint 

(Jy km/s) 


(km/s) 


(km/s) 


(6,6) 


25056.03 


408 


2010/10/11 


K 


0.93 X 0.76 


3.3 


2.54 


30.78 


5.31 


4.85 


(7,7) 


25715.18 


539 


2010/10/12 


K 


0.91 X 0.72 


3.3 


1.62 


16.25 


5.14 


4.01 


(8,8) 


26518.91 


687 


2010/12/17 


Ka 


0.91 X 0.79 


18. 


1.30 


10.41 


4.21 


3.21 


(9,9) 


27477.94 


853 


2010/12/21 


Ka 


0.82 X 0.66 


3.1 


0.88 


6.89 


5.37 


3.15 


(10,10) 


28604.74 


1036 


2010/12/29 


Ka 


0.84 X 0.65 


3.0 


0.61 


3.89 


5.68 


2.54 


(11,11) 


29914.49 


1238 


2011/01/10 


Ka 


0.77 X 0.65 


2.9 


0.34 


1.98 


5.90 


2.31 


(12,12) 


31424.94 


1456 


2011/02/05 


Ka 


0.77 X 0.63 


3.5 


0.17 


0.82 


5.93 


1.91 



^fote. — 

^— Transitions include ortho-NHs {K — 3n) and para-NHs {K ^ 3n). 
^—Energy above the ground from the JPL database. 

RMS noise in a ~0.2 km s~^ channel (no velocity-smoothing). 

Fpcak, Fint! and 17 are estimated from single-Gaussian fits to the spectral profiles in Figure [Tl 
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Figure 2. Total intensity images of NH3 inversion transitions from (6,6) to (12,12) as observed toward Orion BN/KL witfi the EVLA. 
The images were integrated over the velocity range -14 km s~^ to 24 km s~^ with a flux cutoff of 5 mjy beam~^ (~ Scr) for all transitions. 
Contours are 1, 20, 50,.., 550 (by 50) mJy beam~^ km . The beam (0f'85x 0"7) is shown in the right-lower corner. In every panel, the 
crosses mark the position of 4 YSOs in the region: sources I, SMAl, n, and IRc7 (left to right). 



is not surprising considering that these transitions arise 
from levels >850 K above the ground, and are thus ex- 
pected to arise from compact regions close to the exciting 
source. 

We produced total intensity images of various NH3 
inversion transitions, integrated over the whole band 
(Fig. [2]). In the image we included also the positions 
of the four YSOs Source I, n, SMAl, and IRc7. The 
low-excitation NH3 emission (6,6; 7,7) is extended over 



- 15" X 15" (0.03 X 0.03 pc) and shows the typical 
Hot Core "heart-shaped" structure previously observed 
in var ious mo lecular lines: C2 H5CN, CH,30H , OCS 
(jFriedel fc Snvd cr 20081 CH^C N (IZaoata et al.llMTbl) . 
and NH3 (4,4) (Wilson et alj I2000D . The northwest 
(NW) lobe has the weakest integrated emission and it 
is not detected in the higher- J transitions, whose emis- 
sion originates from the northeast lobe of the heart- 
shaped structure. This confirms earlier results from 
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iWilson et all (|2000l ) based on a comparison between (4,4) 
and (10,9). Interestingly, the highly-excited NH3 line 
emission shows an "arc-like" ridge, oriented northeast- 
southwest (NE-SW), <3000 AU long and <1000 AU 
across, curving around Source I. Source I lies at the NW 
edge, indicating a lower column density towards its po- 
sition. 

West of Source I, we find high-J emission associated 
with IRc7. The structure of this emission in the tran- 
sitions from (10,10) to (12,12) is bipolar and elongated 
north-south (N-S), with size -1000 AU. 

3.1. Temperature/density analysis 

We estimated rotational temperature, Tj-ot, and col- 
umn density, Ncoi, from rotational energy diagrams, us- 
ing six observed transitions. We assumed that the inver- 
sion transitions are thermalized and optically thin. The 
beam-ave raged Ncoi from an opt ically thin transition is 
given by (jFriedel fc Snvdeii 12001 : 



^g-£„/T.„* ^ 2-04iy 



X 10^°cm"2 



(1) 



where Q is the partition function, is the upper state 
energy of the transition (K), W is the integrated line 
intensity (Jy beam"'^ km s~^), 9m and 9m are the 
FWHM Gaussian synthesized beam dimensions (arcsec- 
onds), Sji^ is the product of the line strength and the 
square of the molecular dipole moment (Debye^), and v 
is the transition frequency (GHz). For H2 densities in the 
Hot Core (- 10^ cm'^; [Gen zcl fc S tutzkil [l989l) we as- 
sume Ti.ot=Tkin (jWilson et al..,2000^ . so in the following 
we simply refer to Tkin- 

We performed separate temperature analyses for ortho 
and para species, enabling empirical determination of the 
NH3 ortho-para ratio in the Hot Core: 0.9-1.6 (±0.05 - 
0.1) with an average over all pixels of 1.1. The highest 
values are located in the NW-side of the ridge which faces 
Source I. 

We determined the distribution of Tkin and Ncoi on 
a position-by-position basis. We measured Tkin'^lBO- 
490 K (AT < 8 K) and Ncoi-- (1 -4) x 10^^ cm-2(AA^ < 
10^^ cm~^), averaged on scales ^l"(Fig.[3]). The highest 
values of Tkin and Ncoilie in a NE-SW elongated ridge 
offset (to the southeast) from Source I. At Source I's po- 
sition, we measure Tkin~260 K and Ncoi~ 7 x 10^^ cm~^. 
Temperature and density towards Source I are estimated 
using only ortho transitions because the signal-to-noise 
ratios in para data were insufficient. Consequently, total 
c olumn density is pro bably 1.5-2 x greater. 

IWilson et all ()1993D reported 400 ± 40 K using the in- 
version line s (10,10) to ( 1 4,14) observed at 30" resolu- 
tion, while IWilson et al.l ()20 00) reported values in the 
range 130-170 K from 1" images of (4,4) and (10,9). In 
the paper presented here, high-excitation transitions ob- 
served at sub-arcsecond resolution enable us to map the 
gas density and temperature distribution with unprece- 
dented precision. This has revealed hotter, denser ma- 
terial than previously derived from NH3 measurements. 
The tempe rature inferred from high-J lines of CH3CN is 
200-600 K (|Wang et al.ll2010D . consistent with our find- 
ings from NH3. 

Our ana lysis assumes optically thin transitions 
(jWilson et al..i2000i ). Hyperfine-structure satellite lines 



are outside of the bandwidths employed, hence direct es- 
timati on of optical depth is impossible. iHermsen et al.l 
(|1988D reported optical depths >1 for the (6,6) and (7,7) 
transitions. Non-negligible optical depths in the lowest 
available transitions (6,6 for ortho; 7,7 for para) might 
systematically boost the derived temperatures. We ar- 
gue otherwise that this is not the case. Para-transitions 
yield higher temperatures than ortho-transitions, incon- 
sistent with expectation of lower optical depth for (7,7) 
than (6,6). Moreover, rotation diagrams for individual 
pixels do not diverge from a linear trend for the lower 
(J,K) states. 

3.2. Velocity field 

The first moment maps of (6,6) and (7,7) in Figure |4] 
show that redshifted gas (7 to 14 km s~^) is located on 
the NW part of the heart-shaped structure, blueshifted 
gas (—10 to km s~^) is concentrated towards the cen- 
tral part, while the NE lobe is around systemic veloc- 
ity (0 to 7 km s~^). The NH3 emission shows a veloc- 
ity gradient southward from Source I towards the Hot 
Core center, where the largest values of velocity disper- 
sion are also observed. Unless we assume the presence 
of an additional embedded source in that location, this 
may indicate shock propagation in the Hot Core (see 
Sect. [4]). Interestingly, we identified a N-S velocity gra- 
dient (Au '--^5 km s~^) in the bipolar structure traced 
by high-J NH3 emission in IRc7, where line- widths also 
appears enhanced. The N-S bipolar morphology, the ve- 
locity gradient in the elongation axis, and the line-width 
enlargement provide strong evidence of a low-velocity 
outflow associated with IRc7. 

4. DISCUSSION & CONCLUSIONS 

Our multi-transition measurements of NH3 emission 
enable the estimation of temperature and density of (hot) 
molecular gas at high-spatial resolution in the vicinity 
of Source I and the Hot Core. This complements the 
temperature distribution as a function of radius close 
to Source I inferre d from radiative t ransfer modelling 
of SiO masing gas (iGoddi et al.ll2009l ). Maser action in 
the disk/wind from Source I implies temperatures 1000- 
2000 K inside ~100 AU and temperatures of 400-1000 K 
in the outflow at radii 100 AU to 1000 AU. A dust color 
temperatu re of f-700 K is consistent with this flnding 
(Okumura et al.l l 2011| ). On the other hand, an NH3 tem- 
perature of ~260 K toward Source I implies that hot NH3 
emission is not directly (radiatively-) excited by Source I. 
One explanation of this apparent discrepancy could be 
that a high ultraviolet radiation fleld in the vicinity of 
Source I may dissociate fragile molecules like NH3 (dis- 
sociation energy 4 eV), while enhancing gas-phase abun- 
dance of SiO through grain sputtering. This would also 
explain the offset between Source I and the peak of high- 
J NH3 emission. 

What is exciting the high-energy transitions of NH3? 
In fact, all known or suspected (massive) YSOs in the re- 
gion are offset from peaks in hot NH3 emission (Figure[2]), 
and there is no evidence of any other embedded ob- 
jects heating the Hot Core in the radio (iMenten fc ReidI 
I1995D. millimeter (iFriedel fc Snvde^ l2008D. and infrared 
(|GreenhiU et al.l l2004bD . Additionally, the color tem- 
perature distribution from mid-IR measurements shows 
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Figure 3. NH3 rotational temperature (Trot) and column density (Ncoi) maps of Orion BN/KL, obtained from ortho (colors) and para {contours 
in top two panels) inversion transitions. In all panels, the crosses mark the positions of sources I, SMAl, n, and IRc7 (left to right), (a) Map of 
Trot averaged over LSR velocities —14 km s~^ to 24 km s"^. Contour levels are 210 K to 480 K in increments of 30 K. Outside the region shown, 
the signal-to-noise ratio was insufficient to calculate Trot (AT >8 K) reliably, (b) Map of Ncoi obtained using equation 1. Contours are from 0.8 
to 2 xlO"*""^ cm~^ in increments of lO"*"^ cm^"^ {A.N < 10^^ cm""^). Note that the densest and hottest NH3 are found in a northeast-southwest 
elongated ridge to the southeast of Source I. (c) Ortho N^oi (colors) with SiO v—0, J — 2 — 1 total intensity (contours) from fPlambeck et al.| ||2009[) 



overplotted. Conto urs are 0.1, 0.3, 0.5, 0.7, 1. 
jGoddi et aI.|[20TT|) . 



2, 3 Jy beam km s . The arrow indicates the proper motion of Source I in the Orion rest-frame 
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Figure 4. Velocity fields of NH3 inversion transitions from (6,6) to (12,12). Colors indicate Y lsr in km s~^. Color scales are compressed 
with increasing quantum number to clearly show velocity structure. In each panel, the crosses mark the position of sources I, SMAl, n, 
and IRc7 (from left to right). All transitions show a velocity gradient southward from Source I towards the Hot Core center (the less clear 
signature in the upper panels is due to the expanded color scale). 
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a decreasing gradient from Source I to the Hot Core 
(|Okumura et aLllMll) . We propose that the elevated 
gas temperatures measured from NH3 in the Hot Core 
are produced by shock-waves and not by steUar ra- 
diation. Theoretical models suggest that NH3 might 
form through surface chemistry on dust grains and 
then be released in the gas-phase by the passage of 
shocks, either through hydrodynam ic gas-grai n inter- 
actio n (|Flower fc Pineau des Foreti ll994; Flow er et al 



1995 ) or shock- induced IR radiation (jTavlor fc Williams 
1996( ). Indeed, enhanced NH3 has been observed in 



the L1157 low-mass outflow (Tafall a fc Bachillei] IT995h 
as well as towards a dense (quiescent) clump down 
stream from the hi gh -mass protostellar jet H H 80 North 
(jGirart et al.lll998D . lUmemoto erall (|19990 report also 
ortho-para enhancement in the L1157 outflow. Several 
pieces of evidence corroborate the shock hypothesis in 
the Orion Hot Core: the offset of hot NH3 emission 
peaks from known (proto)stellar sources, the high Tkm, 
the departure of the ortho- para ratio from the st atisti- 
cal equilibrium value fl.O: lUmemoto et al.lll999l ). and 
the high abundances of shock tracer molecules (SiO, SO, 
SO2) and common molecular con stituents of interstellar 
ices (C2H5CN, CH3OH, CH.CN: Iro^ Dishoeck fc BlaO 
[19981) . 

What generates the shocks? Recent studies report 
strong evidence that a dynamical interaction occurred 
^^500 years ago betweeii Source I and the hi gh-mass YSO 
BN (Gom ez et al1[200l ICoddi et al.l[20Tll ). The stellar 
interaction resulted in high stellar proper motions for 
Source I and BN and possibly pro duced the fast bul- 
let outflow traced by CO and H2 (jZapata et all 120091: 
iBallv et al.ll20TTI: iGoddi et al.ll201ll ). In this framework, 
it has recently been proposed that the Hot Core might 
have originated from the impact of a shock-wave onto 
a pre-existing dense core, driven by either the expan- 
sion of a dense bubble of material previousl y associated 
with BN an d Source I (iZapata et al.ll2011a|) and/or the 
fast bullets (jZapata et al.l l2011b[ ). Data in Figure |3]:, 
however, suggest another possibility. The relationship 
between the high column density gas as traced by NH3, 
the Source I proper moti on, and the Source I o utflow as 
traced by SiO emission (jPlambeck et al.ll2009D together 
suggest that the Source I outflow drives compression and 
shocks in the surrounding medium. The relative narrow- 
ness of the SW outflow in close proximity of the column 
density peak is suggestive of ram pressure effects caused 
by the outflow and the YSO proper motion. The more 
open NE lobe, its apparent symmetry about the flow 
axis, and its extension beyond the northern NH3 clump 
may indicate that the outflow is in front of or behind the 
dense gas, compressing it more along the line of sight. 
Part of the mechanical energy of the outflow could be 
dissipated and contribute to gas heating. This is also 
consistent with the ortho-para ratio being larger in the 
NW-side of the ridge facing Source I. 

Detailed modeling, including shock-physics as well as 
dust/gas-phase chemistry, will be required to discrim- 
inate between different (shock) scenarios proposed for 
the Hot Core excitation. However, the shock velocities 
implied by different scenarios and inferred NH3 chem- 
istry are suggestive. iFlower fc Pineau des Foreta (|1994D 
showed that shocks with velocities of 10-50 km s"-'^ (C- 



type) efficiently sublimate ice mantles, ejecting volatile 
molecules intact and enabling efficient abundance en- 
hancement. Shocks with velocities <10 km s~^ would 
have just a diminished impact on grain-mantle evap- 
oration and NH3 abundance enhancement, which 
weighs against the hypothesized slowly-expanding bub- 
ble (7 km s~^j3 as the main mechanism for the Hot 
Core excitation. In contrast, fast bullets (> 100 km s""'^) 
would drive strong J-type shocks, causing grain sputter- 
ing and molecular dissociation. Reformation of NH3 by 
gas-phase chemistry is possible but less likely on the short 
dynamical time of the region (^500 years). Nonethe- 
less, existing models have not investigated chemical evo- 
lution behind J-type shocks. In termediate-velocit y bul- 
lets traced by H2 (30-70 km s~^: IBallv et"all[20n could 
in principle drive C-type shocks and release NH3 as well 
as excite H2 emission, but the poor correlation betw een 
NH3 and mid-IR/H2 emission l|Shuping et al.ll2004l ) ar- 
gues against this possibility. We note that the impact of 
outflow from Source I (20 km s~^) in combination with 
proper motion of 12 km s~^could efficiently h eat the Hot 
Core with limited molecular dissociation. iFavre et al.l 
(|2011h reported an anticorrelation of methyl- formate and 
H2 in the Hot Core and proposed a similar mechanism for 
its excitation. Extinction could however be responsible 
for the observed anticorrelation. 

In summary, we suggest that local mechanical interac- 
tion of Source I with the Orion Hot Core is more likely to 
be responsible for enhanced NH3 abundance and heating 
of gas in the Core than the consequence of the hypoth- 
esized explosive event tied to the protostellar dynamical 
interaction. 
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